A Particle Model of Rolling Grain Ripples Under Waves 
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A simple model is presented for the formation of rolling gram ripples on a flat sand bed by the 
oscillatory flow generated by a surface wave. An equation of motion is derived for the individual 
ripples, seen as "particles", on the otherwise flat bed. The model account for the initial apperance 
of the ripples, the subsequent coarsening of the ripples and the final equilibrium state. The model is 
related to physical parameters of the problem, and an analytical approximation for the equilibrium 
spacing of the ripples is developed. It is found that the spacing between the ripples scale with the 
square-root of the non-dimensional shear stress (the Shields parameter) on a fiat bed. The results 
of the model are compared with measurements, and reasonable agreement between the model and 
the measurements is demonstrated. 
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In the coastal zone where the water is relatively shal- 
low, an ubiquitous phenomenon is the formation of rip- 
ples in the sand. These ripples have been described in 
the seminal work of Bagnold 0] , who called them vortex 
ripples. The vortex ripples have been studied recently 
from a pattern forming point of view . Bagnold also 
described another kind of ripples which were created only 
from an initially flat bed, as a transient phenomena be- 
fore the creation of vortex ripples. These ripples were 
created by the rolling back and forth of individual grains 
on the flat bed, and were thus called rolling grain ripples. 
It is this latter class of ripples which is the topic of this 
article. 

The most well-known type of ripple created by oscil- 
latory motion is the vortex ripple, called so because of 
the strong vortices created by the sharp crest of the rip- 
ples. These are the ripples one encounter when swimming 
along a sandy beach. The shape of the ripples is approxi- 
mately triangular, with sides being at the angle of repose 
of the sand. The length of the ripples scale with the am- 
plitude of the oscillatory motion of the water near bed, 
a. The dynamics of the vortex ripples and the creation of 
a stable equilibrium pattern has been recently described 

by §• 

In contrast to the vortex ripples, the rolling grain rip- 
ples are rarely (if ever) found in the field, and have only 
been studied in controlled laboratory experiments, where 
it is possible to have a completely flat bed as initial con- 
ditions. In this article the term "rolling grain ripple" 
refer to the ridges with triangular cross section occur- 
ring on an otherwise flat bed (Fig. [I] and Fig. ^) which 
is in agreement with Bagnolds definition of the rolling 
grain ripples (see also Plate 3 in jjj]) . The length of the 
triangular cross section is small (typically smaller than 
1 cm), and much smaller than the spacing between the 



ridges, which again is typically smaller than the spacing 
between the vortex ripples. As the crest of the rolling 
grain ripples is rather sharp, a separation zone can be 
induced in the lee side of the ripple (a separation zone 
is an area where the flow is in the opposite direction of 
the main flow) . This has also been found from numerical 
simulations of the flow over triangular ridges Q . 

Bagnold assumed that there were no separation zones 
behind the rolling grain ripples, and this has lead to the 
adoption of the term rolling grain ripples for ripples with- 
out separation. It seems as if another type of ripple can 
emerge from a flat bed. These ripples initially have a si- 
nusoidal shape, and thus no flat bed between them. The 
ripples grow and become sharper until they are steep 
enough to induce separation, whereafter they grow to 
become vortex ripples. From the authors own experience 
H and from |7j], it seems as if which if these two kinds 
of ripples appear from a flat bed depend on the initial 
preparation of the bed. If the grains are small and the 
bed well packed, the rolling grain ripples as described by 
Bagnold appear ("Bagnold type rolling grain ripples"). 
If, however, the grains are large or the grains are care- 
fully prepared, so as to not be packed, i.e., by letting the 
grains settle trough the water, a mode with sinusoidal rip- 
ples can be observed 0. These "sinusoidal rolling grain 
ripples" have been described by a hydrodynamic instabil- 
ity of the wave boundary layer. The idea was originally 
formulated by Sleath ||, but was developed in great de- 
tail in a series of papers by Blondeaux, Foti and Vittori 
|^| |l3| . This kind of ripples is not described by the model 
which will be developed here. 

The model to be presented is a "particle model" in the 
sense that it view each ripple as a "particle" which in- 
teract with its neighbouring particles. A heuristic equa- 
tion of motion is written for each particle, taking into 
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account the movement of the particle and the presence 
of the neighbouring particles. When two particles col- 
lide they merge and form a new, larger particle. This 
continues until a steady state is achieved. 

After a more detailed description of the rolling grain 
in section [| the heuristic model is developed in section 
[n|. It turns out, fortunately, that the parameters enter- 
ing this mod el can be co nnect ed to measurable quantities 
(section II A). In section II B the model is solved numeri- 
cally and analytically, and comparison with experimental 
mea surements of rolling grain ripples are done in section 



II C . Conclusions can be found in section IV 



I. THE ROLLING GRAIN RIPPLES 

An example of a flat bed with rolling grain ripples co- 
existing with vortex ripples is seen in Fig. |l|. In the mid- 
dle of the picture the ripples have not yet formed and the 
bed is still flat, while on the top vortex ripples are seen to 
invade into the flat bed. The vortex ripples are typically 
nucleated from the boundaries or from a pertubation in 
the bed. In the lower part of the picture the rolling grains 
ripples have formed on the flat bed, and are seen as the 
small bands of loose grains on top of the flat bed. 

The flow over the bed created by the surface wave is 
oscillating back and forth in a harmonic fashion. This 
flow creates a shear stress on the bed r(t) which in non- 
dimensional form reads: 



9{t) 



r(t) 



p(s - l)gd 



(1) 



where p is the density of water, s is the relative density 
of the sand (for quartz sand in water s = 2.65), g is the 
gravity and d is the mean diameter of the grains. 9 is usu- 
ally called the Shields parameter [l4| . When the shear 
stress exceeds a critical value 8 Cl the grains start to move. 
For a turbulent boundary layer the value of the critical 
Shields parameter is 9 C sa 0.06 |L5|. The grains which 
have become loosened from the bed, start to move back 
and forth on the flat bed, and after a while the grains 
come to rest in parallel bands. In the lee side of each 
band the bed is shielded from the full force of the flow, 
creating a "shadow zone" where the grains move slower 
than in the upstream side of the bands. 

Due to this shadow zone, more grains end up in the 
bands than leave the bands, and they grow until they 
form small ridges, the rolling grain ripples. When the 
rolling grain ripples are fully developed, no grains will be 
pulled loose from the bed in the space between them, and 
they are stable. However, in reality the rolling grain rip- 
ples are dominated by invading vortex ripples, which is 
the main reason why they are rarely observed in Nature. 



II. A SIMPLE MODEL 

The above scenario can be formulated mathematically 
by writing an equation of motion for each grain/particle. 
In the following an equation of motion for the particles 
are developed. In the beginning the particles represents 
the grains, but as the single grains quickly merge, the 
particles most of the time represents a ripple. First the 
velocity of each particle is found, assuming that the par- 
ticle is alone of the flat bed, and then the influence of the 
shadow zones from neighboring particles are taken into 
account. 

Consider N particles rolling on top of a rough, solid 
surface. Each particle is characterized by its position Xi 
and its height hi (see Fig. ||). As the ripples are trian- 
gular, the area of each particle Ai and their heights are 
related as: 



hi = y Ai tan i 



(2) 



where <f> is the angle of repose of the sand (33°). 

The ripple moves back and forth slower than a sin- 
gle grain, according to the "1/height" law. This law is 
well known in the study of dunes in the desert |l6[] or sub- 
aqueous dunes 17 1, and can be illustrated by a simple ge- 
ometrical argument. Suppose there is a flux of sand over 
the crest of a ripple or a dune q cres t (Fig. ||). To make 
the ripple move a distance 5x an amount of sand hSx is 
needed. As the sediment flux is amount of sand per unit 
time, the velocity of the ripple is u r i pp i e — q cres t/h oc 1/h. 
If the height of the initial particles (the single grains) is 
assumed to be equal to the grain diameter d, the velocity 
of the particles can be related to that of the single grains 
as: 



Ui = — U g sin(wt) , 

tli 



(3) 



where U g is the velocity amplitude of the motion of a sin- 
gle grain and ui is the angular frequency of the oscillatory 
motion. 

The flow moves back and forth and makes the particles 
move accordingly on top of the bed. In the wake of each 
particle/ripple there is a shadow zone (Fig. 0), which is 
the area behind the particle where the absolute value of 
the shear stress is smaller than it would be on a flat bed. 
The length of the shadow zone is therefore larger than 
the length of the separation bubble formed by the parti- 
cle (note that the shadow zone would be present even in 
the absence of separation). If the shadow zone is much 
smaller than the amplitude of water motion, the flow in 
the lee side of the ripple can be assumed to have sufficient 
time to become fully developed. The fully developed flow 
over a triangle is similar to that past a backward facing 
step in steady flow, which has been extensively studied 
(see e.g. Tjerry fig]). In that case the relevant quanti- 
ties, i.e., the length of the separation bubble, the length 
of the wake etc., scale with the height of the step. As a 
first assumption, the shadow zone is therefore assumed 
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to have a length which is proportional to the height of 
the particle: a s hi. If a particle enters the shadow zone of 
another particle, it is slowed down according to the dis- 
tance between the particles. This means that the actual 
velocity of a particle is Uif(Ax) where Ax is the distance 
between the grain and the nearest neighbor upstream, Ui 
is the velocity of the particle outside the shadow zone 
and / is a function determining the nature of the slowing 
of the motion of the particle. A simple linear function is 
used, as shown in Fig. ||. The exact form of the function 
/ is not crucial, as will be evident later, the important 
parameter being the extent of the shadow zone as deter- 
mined by a s . 

It is now possible to write the equations of motion for 
the particles as a system of coupled ODEs: 

^ _ d_ u u\f( Xi Z Xi ~ 1 u ^ \ ff Xi ~ Xi+1 u ^ \ 
hi 9 \ a s hi-i \u(t)\J \ a s h i+ i |tt(i)|/ 

S v /V v ' 

positive half period negative half period 

(4) 

for i = 1,...,N, where u g (t) = U g sm(ujt). The mo- 
tion of a particle is thus made up of three parts (i) the 
motion of the single undisturbed particle, (ii) the effect 
of the shadow from the particle to the left (i— 1), which 
might affect particle i in the first half period and (Hi) the 
effect of the shadow of the particle to the right (i + 1) in 
the second half period. 

When lengths are scaled by the diameter of the grains 
and time by the frequency lj, it is possible to identify the 
three relevant dimensionless parameters of the model: 

a s , the length of the shadow zone divided by the height. 

a g /d, the amplitude of the motion of a single grain, di- 
vided by the grain diameter (a g — ioU g ). 



1. The velocity of the grains 



The velocity of the grain can be found by considering 
the force balance on a single grain lying on the flat bed. 
The grain is subject to a drag force proportional to the 
square of the relative flow velocity u r = u n b — u g where 
u n b is the velocity near the bed and u g is the velocity of 
the grain: 

Fd = ^C D pA\u r \u r , (5) 

where A is the area of the grain and Co is a drag coeffi- 
cient. The numerical signs is used to obtain the right sign 
of the force. The velocity profile in the vicinity of the bed 
is supposed to be logarithmic. As shown in po) this is 
a reasonable assumption except when the flow reverses. 
However, during reversal the velocities are anyway small 
and the accuracy is of minor importance. The logarith- 
mic profile over a rough bed can be written as (eg. Jul): 

«(y) = ^ln(^), (6) 
K \k N J 

where k = 0.41 is the von Karman constant, Uf = \frjp~ 
is the friction velocity and is the Nikuradse roughness 
length. It is then possible to find the near bed velocity 
as the velocity at yml/2d: 

Unb = frf (7) 

where the constant £ can be determined from Eq. (^|) by 
assuming fcjv = d. Opposing the drag on the grain is the 
friction of the bed: 



Xi/d, the initial distance between the grains; Aj = L/N, 
where L is the length of the domain and N is the 
initial number of grains. 



A. Relation to physical quantities 



F f - - f iW 



(8) 



where fi is a friction coefficient and W = pg(s — 1)g? 3 7t/6 
is the immersed weight of the grain. By making a bal- 
ance of forces, Fd + Ff = 0, the velocity of the grain can 
be found: 



Even though the model seem quite heuristic, the pa- 
rameters entering the model, a g /d, Xi/d and a s can be 
related to physical parameters describing the flow and 
the properties of the grains. The line of arguments pre- 
sented here closely follows those used to derive the flux 
of sand on a flat bed (the bed load), as can be found in 
text books, i.e., |15| , or in [|9|. First the velocity of a 
single grain will be derived, from which a g /d can be in- 
ferred. Thereafter the initial number of grains in motion 
is found, from which follows Xi/d. Finally the length of 
the shadow zone a, is discussed. 



where 



£u f 1 



4/i 



is the critical Shields parameter. 
0.65 [EH. 



(9) 



(10) 



Usually fi = tan< 
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2. The initial spacing of grains 

Now that the velocity of the grains has been calculated, 
there still remains to determine the number of grains per 
area n in motion. To this end a small volume of moving 
sand at the top the flat bed is considered. The balance 
of the forces acting on this volume is written as: 



n 



(11) 



The interpretation of the terms is as follows: The param- 
eter Tb is the shear stress on the top of the bed load layer. 
It is assumed that this is equal to the shear stress on a 
fixed flat bed. tg is the stress arising from the inter- 
granular collisions, giving rise to "grain-stresses" Jl9| , 
modeled as: tq = nfiW. It is assumed that the inter- 
granular stress absorbs all the stress, except the criti- 
cal stress t c ; this is the so-called "Bagnold hypothesis" 
|L9j . Making Eq. (|ll]) non-dimensional by dividing with 
p(s — l)gd, the number of grains in motion is found as: 



G 



(12) 



If 6 < 6 C , then no grains are in motion and n = 0. n can 
also be viewed as the initial density of grains, and by as- 
suming a square packing of the grains the initial distance 
between the grains becomes: 



A* 
d 



1 



i/nd 



7T ft 



6(e - e c ) ' 



(13) 
(14) 



3. The length of the shadow zone 

The last parameter a s , which characterize the length 
of the shadow zone, is a bit harder to estimate accu- 
rately. By exploring the analogy with the backward fac- 
ing step which was suggested in section H, we can get 
some bounds for a s . In the backward facing step there is 
a zone with flow separation which extends approximately 
6 step heights from the step, see Fig. H. After approxi- 
mately 16 step heights there is a point where the shear 
stress has a small maximum. It follow that the length 
of the shadow zone should be longer than the separation 
zone, but shorter that the point of the maximum in the 
shear stress, i.e., 6 < a s < 16. 



B. Numerical and analytical solutions of the model 

In the following section the behavior of the model is 
examined. To study the detailed behavior the set of cou- 
pled ordinary differential equations (|J) is integrated nu- 
merically. It is seen that the model reaches a steady state 



and an analytical expression for the spacing between the 
ripples in the steady state is developed. 

The numerical simulations in this section are based on 
a simple example with a g /d — 35 and Xi/d = 3.23 and 
a s is set to 10. As initial condition all particles have an 
area of 1.0 ± 10 %, to add some perturbation. The initial 
number of particles N in this example is 800. 

In the first few periods a lot of grains is colliding and 
merging (Fig. |J). As the ripples are formed and grow 
bigger the evolution slows down, until a steady state is 
reached (Fig. |). The spacing between the ripples in the 
steady state show some scatter around the average value, 
which is also observed in experiments. The variation in 
the average spacing of the steady state X eq for realiza- 
tions with different initial random seed, turned out to be 
on the order of l/N eq where N eq is the final number of 
ripples. 

To study the behavior of the average spacing of the 
ripples X eq , a number of simulations were made where 
the parameters were varied one at a time. Each run was 
started from the initial disordered state. 

Changing a g /d only results in a minor change in the 
spacing of the ripples (Fig. ||a). The final spacing be- 
tween the ripples does depend on the length of the 
shadow zone a s ; the longer the shadow zone, the longer 
the ripples (Fig. ||b). This can be used to estimate the av- 
erage equilibrium spacing between the ripples. When the 
distance between two ripples is longer than the shadow 
zone of the ripples, they are no longer able to interact. 
This gives: 



X eq > a„h 



eq j 



(15) 



where subscript eq denotes average value at equilibrium. 
However if two ripples have a spacing just barely shorter 
than Eq. (JlJ), they will be able to interact and eventu- 
ally they will merge. One can therefore expect to find 
spacings up to X eq <2a s h eq . Assuming that the average 
length is in between the two bounds one get: 



A, 



1 < 7 < 2 



(16) 



where 7 can be found by comparing the results from the 
full simulations with Eq. (|l6|). The height of the rip- 
ples at equilibrium can be found by splitting the initial 
number of particles evenly onto the equilibrium ripples. 
Then the average area at equilibrium is A eq = X eq /Xid 2 
and from Eq. (0) follows that the height is h eq /d = 
•\/tan</> X eq /Xi, which gives an average equilibrium spac- 
ing: 



Ae 



2 ajd tan < 



7 



2 2 
a s 7 



^^Ve~e- C . (17) 



The equilibrium spacing is therefore found to be propor- 
tional to \/0 — 9 c with the constant of proportionality be- 
ing made up of a s , 7 and various geometrical factors. All 
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the quantities related to the dynamical evolution of the 
ripples, i.e., the velocity of the ripples, the shape of the 
function /(Ax) etc. does not enter into the expression. 



C. Comparison with experiments 

The only parameter that has not been accurately de- 
termined is a s . The value of this parameter can be esti- 
mated by comparison with measurements. 

In 1976, Sleath made a series of experiments, measur- 
ing spacing between rolling grain ripples |8| . The ripples 
were formed on a flat tray oscillating in still water using 
sand of two different grain sizes: 0.4 mm and 1.14 mm. 
To compare with the experiments the value of 9 need 
to be calculated. 9 reflect the number of grains in mo- 
tion, and it seems natural to use the maximum value dur- 
ing the wave period, 9 max . To estimate 9 max the shear 
stress on the bed has to be estimated. The maximum 
shear stress on the bed, T max during a wave period can 
be found using the concept of a constant friction factor 
U ( @): 



_1 . . 

7~max — ^PJw^ 



2 

max "> 



(18) 



with u max being the maximum near-bed velocity. The 
friction factor can be estimated using the empirical rela- 
tion £y|: 



'» = °- 04 { M 



-0.25 



(19) 



where a 



r u> and k 



N 



d. 



The range of Shields parameters in the experiments 
was found in this way to range from the critical Shields 
parameter to 9 — 0.42. For the high Shields parame- 
ters the rolling grain ripples were reported to be very 
unstable, and to quickly develop into vortex ripples. In 
these cases, the measured ripple spacing then reflects the 
spacing between the rolling grain ripples before they de- 
veloped into vortex ripples M . 

In Fig. ^| the experimental results are compared with 
runs of the model using a s = 15.0 (the reason for this 
particular value will be clear shortly) and N = 10000. By 
fitting all the runs to Eq. d|) it was found that 7 = 1.40. 
The results using Eq. ( p"7j ) and 7 = 1.40 are shown with 
a line. 



First of all, it is seen that Eq. (17) predicts the results 
from the full model Eq. (|J) well. The correspondence be- 
tween the model and the experiments is reasonable, but 
there are some systematic discrepancies, which will be 
discussed. 

There are some few points with small ripple spacing for 
which the model does not fit the measurements. These 
measurements has a Shields parameter very near the crit- 
ical (i.e., just around the onset of grain motion), which 
implies some additional complications. The grains used 



in the experiment were not of a uniform size; rather they 
were part of a distribution of grain sizes, and the grain 
size reported is then the median of the distribution, d^. 
The Shields parameter is calculated using the median 
of the distribution, but actually one could calculate a 
Shields parameter for different fractions of the distribu- 
tion, thus creating a 9\o, a #50 etc.. When #50 is smaller 
than then the critical Shields parameter, 9\q might still 
be higher than the critical Shields parameter. This im- 
plies that grains with a diameter smaller than d$o will 
be in motion, while the larger grains will stay in the 
bed. As only d$o is used in the calculation of the equilib- 
rium ripple spacing, the distance between the grains A; 
will be overestimated near the critical Shields parameter, 
where the effect of the poly-dispersity is expected to be 
strongest. An overestimation of A; will lead to an under- 
prediction of the ripple length, which is exactly what is 
seen in Fig. |[ 

There are also three points from the experiments taken 
at very large Shields parameters, which seems to be a bit 
outside the prediction of the model. As already men- 
tioned, these points are probably doubtful because of the 
very fast growth of vortex ripples. It is therefore rea- 
sonable to assume that vortex ripples invaded the rolling 
grain ripples before these had time to reach their full 
length. 

To find a reasonable value of a s Eq. ([l?]) was fitted 
to the experimental points. To avoid the points which 
might be of doubtful quality, as discussed above, only 
the points in the range 0.075 < 9 < 0.3 were used. This 
gave the value of a s = 15.0, in agr eement with the quali- 
tative arguments in section [I A 3. 



III. DISCUSSION OF THE RESULTS 



From the comparison of the model with measurements 
it seems as if the model confidently reproduces the ex- 
periments. 

In the model the number of grains in motion is con- 
stant (even though the number of particles changes). In 
an experimental situation, however, new grains might be 
lifted from the bed and added to the initial number of 
grains in motion. As the part of the flat bed between 
the ripples is covered by the shadow zones of the parti- 
cles, these stretches will be shielded from the full force 
of the flow, and only very slowly will new grains be loos- 
ened here. This small addition of new grains will result 
in a slow growth of the rolling grain ripples, such that 
eventually grow into vortex ripples. This slow growth is 
very well illustrated by recent measurements Q , but not 
covered by the present model. 
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IV. CONCLUSION 

In conclusion, a model has been created which explain 
the creation and the equilibrium state of rolling grain 
ripples of the type described by Bagnold. The final dis- 
tance between the ripples is proportional to y/6 — 6 C . The 
model has been compared with measurements with rea- 
sonable agreement. 
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FIG. 1. An example of a flat bed viewed from above, with rolling grain ripples (bottom) and vortex ripples invading from 
the top. The experimental setup is a wave tank 60 cm wide and approximately 50 cm deep, containing sand with a median 
diameter of 0.2 mm. 




FIG. 2. A sketch of a ripple on the flat bed with related quantities. 
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FIG. 3. Geometrical illustration of the quantities involved in the derivation of the one-over-h law. 
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FIG. 4. An illustration of the shadow zone of one ripple, in the part of the wave period where the flow is from the left to 
the right. Below is seen the function /(Ax) used to describe how a particle is slowed down when it enters the shadow zone of 
another particle. 



6h 16h 

FIG. 5. A sketch of the shear stress on the bed during a steady flow over a backward facing step. 
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FIG. 6. Zoom of the movement of the particle in the first two wave periods. a g /d = 35, A/ = 3.23 and a 3 = 10.0. 
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FIG. 7. The development of the particles until steady state is reached. The lines show the positions of the particles at the 
end of each period. 
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FIG. 8. The spacing between the rolling grain ripples at equilibrium as a function of the initial density of grains. The basic 
example is shown with the circles: a g /d = 35 and a s — 10.0. In a) the value of a g /d is varied, for pluses a g /d = 20, crosses 
a g /d — 50. In b) a s is varied; pluses: a s — 7, crosses: a s — 13. 
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